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Abstract 

The reaction rates of the beta processes for aU particles at arbitrary degeneracy are derived, and 
an analytic steady state equiUbrium condition = /ip + 2/Xe which results from the equality of 
electron and positron capture rates in the hot electron-positron plasma with nucleons is also found, 
if the matter is transparent to neutrinos. This simple analytic formula is valid only if electrons 
are nondegenerate or mildly degenerate, which is generally satisfied in the hot electron-positron 
plasma. Therefore, it can be used to efficiently determine the steady state of the hot matter with 
plenty of positrons. Based on this analytic condition, given the baryon number density and the 
temperature, if the nucleons are nondegenerate, only one algebraic equation for determining the 
electron fraction is obtained, which shows the great advantage of the analytic equilibrium condition. 
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I. INTRODUCTION 

Gamma Ray Bursts (GRBs) ^ and core collapse supernovae(SNe) Q are two of the most 
violent events in our universe. Ironically, their explosion mechanisms are still mysterious. 
Recently, the central engine of GRBs is believed to be related to the hyperaccretion of a 
stellar-mass black hole at extremely high rates from ~ 0.01 to 10 MqS~^ 0,0, In such an 
accretion disk, matter is so dense that photons are trapped. The possible channel for energy 
release is either neutrino emission which is mainly from the electron-positron (e^) capture on 
nucleons and annihilation, or outflows from the disk. Whatever a successful central engine 
is, it ejects a hot fireball which consists of the radiation field (e^/photons) and baryons. The 
ratio of neutrons to protons, or equivalently, the electron fraction, is crucial to the observed 
radiation from GRBs [7|, |8|, its dynamic evolution and the nucleosynthesis in the disk 
or fireball 14,^1. For instance, the inelastic collisions between neutrons and protons produce 
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observable multi-GeV neutrino emission [1|, 110||. In addition, the two component fiuid of 
neutrons and protons significantly changes the fireball interaction with an external medium 
which is su ppo sed to produce the observed electromagnetic radiation from GRBs and their 
afterglows [12]; and the electron fraction strongly affects the equation of state of the 
hyperaccretion disk and the neutrino emissions from it. 

Roughly speaking, SNe are powered by the iron core collapse of their progenitors. Most 
numerical simulations have shown not only the failure of the prompt shock, but the failure 
of its revival by the delayed neutrino emission from the protoneutron star (PNS). The re- 
sult (explosion or not) sensitively depends on the input microphysics, such as the electron 
capture, the neutrino emission, neutrino-matter interactions, the equation of state, rotation, 
magnetic field, general relativity effects, and so on (see ref. ^, and references therein). With- 
out a doubt, weak interactions, especially e"^ captures and neutron decay, play a key role in 
both GRBs and SNe. During the accretion or collapse, these processes exhaust electrons, 
thus decrease the degenerate pressure of electrons. Meanwhile, they produce neutrinos which 
carry the binding energy away to power the explosions. Therefore, electron capture is crucial 
to the formation of the bounce shock of SNe, and the resulting neutrino spectra strongly 
influence the neutrino-matter interactions which are energy dependent and are essential for 
collapsing simulations [iS^ . 

The existence of a hot state with nucleons is the common characteristics of both GRBs and 
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SNe, as well as the PNS, the bounce shock and the early universe jl^. In these systems, the 
beta reactions are the most important physical processes If the system is transparent 
to neutrinos, the steady state is achieved via the following beta reactions [3], 

e~ +p n + Ue, (1) 
+ n — > p + Ue, (2) 
n p + + Ue- (3) 

These beta reaction rates are calculated in the previous studies, usually under one of three 
approximations: the nondegenerate approximation Q|, the degenerate approximation 



20|, and the elastic approximation in which there is no energy transfer to nucleons [21[. In 
ihis paper, applying the structure function formalism developed by Reddy et al. j3](see also 
2^), I derive the reaction rates of the beta processes in the dense subnuclear matter for all 
particles at arbitrary degeneracy. In addition, I find an analytic expression for determining 
the kinetic equilibrium between electron capture and positron capture, which is efficient 
to determine the steady state of the hot matter with plenty of positrons. If the neutrinos 
are partially trapped, neutrino transport should be considered carefully, which is out of the 
scope of this work jl7l|. 



II. REACTION RATES. 

The lowest order Feynman diagrams for reactions are shown in Fig. ^ Because in 

our consideration the energies of leptons are less than a few hundred MeV which is greatly 
less than the rest mass of mesons, thus the interaction Lagrangian from Weinberg- 
Salam theory is reduced to the original Fermi's current-current interaction form. From 
Fermi's golden rule, the reaction rates of the processes (d}-© read (we set h = c = = I) 

A = 2/n 

i=i 

where Pi = {Ei,pl) denotes the four-momentum of particle i {i = z/g/z/g, e'/e"*", n, p). 
Pi = \pi\, and Pi and Pf are the total initial and final momentum, respectively. JF de- 
notes the final-states blocking factor. In reaction (l)-(3), jF=/g/p(l — /„), jF=/g+/„(l — fp) 
, J^=fn{l — — fe) respectively, where fi is the Fermi-Dirac function of particle i. 

In this paper, we just consider the case that the emitted neutrinos can escape freely 
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FIG. 1: The lowest order Feynman diagrams for /? processes, (a) Electron capture by proton: 
e~ + p ^ n + i/g, the energy conservation requires = qo + E^. (b) Positron capture by neutron: 
e"*"+n p+i^e, the energy conservation requires Ef,+ = qo+Ep. (c) Neutron decay: n p+e~ +iJe, 
the energy conservation requires Ef, = —qq — Ep. 

from the system. |Mp in the above equation is the averaged transition rate which de- 
pends on the initial and final states of all participating particles. Screening corrections 
to the electron capture rates in dense astrophysical environments have been investigated 
by many authors In the hot and dense npe^ gas, the Debye radius td is about 

3fm(T/10i°K)i/2(ng/0.01fm-3)-i/2(i;/0.5)-i/2^ where ng = n„(/i„,T) +np(/ip,T) is the 
baryon number density. Thus in our calculations, Coulomb waves are a ppr oximately re- 
placed with plane waves for charged particles. As argued in Reddy et al. j22|, for nonrela- 
tivistic noninteracting baryons {n^ < ^n^nc), where n^^c = 0.16 fm~^ is the empirical nuclear 
equilibrium density, the transition rate averaged over the initial spins becomes a constant, 
i.e. |Mp = G|cos^6'c(l + 3(7^), here Gp ~ 1.436 x 10~^^ erg cm^ is the Fermi weak inter- 
action constant, Oq (sin^^c = 0.231) is the Cabibbo angle, and gp^ = 1.26 is the axial-vector 
coupling constant. For the case of relativistic interacting baryons, the transition rate is 
expressed in terms of the target particle retarded polarization tensor in the relativistic mean 
field theory for interacting baryons, the resulting expression is so complicated (for more 
detail, see 0, |2^). In this paper, the beta reactions at subnuclear density are considered, 
therefore. Taking IMP to be constant is a good approximation. 

n 

Using the structure function formalism developed by Reddy et al. [22], the above integra- 
tions can be simplified into only three dimensional ones, therefore, the rate of the e^-captures 
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Xe-p, Ae+n, and the rate of neutron decay A„ are given by 



Sir 



|M|2 / dE,. 
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'0 Jm^ — E^ 

■\Pe+Pv\ 

X / dqEyEefeSp^n{qo,q)q 

\Pe-Piy\ 



(5) 



X 



Jo Jnie—Eii 

dqEpEe+fe+Sn^p{qo, q)q, 



(6) 



A. 



JO J -oo 



X / dqEpEe{l - fe)Sn-.piqo,q)q, 



(7) 



Pe-Pii| 



where Si^j{qo,q) is the so-called dynamic form factor or structure function which charac- 
terizes the isospin response of the system [22], and qo = Ef — Ei, q = \q\ = \pf — Pi\ denote 
the momentum and energy transfer. 
The expression Si^j{qQ,q) is given by 

rriimfT z + ^- 



Si^j{qo,q) 



Tiq 



exp(— 2;j 



where 



E' 



mi+ 



qo + fJ'i- ftj 
T 

l + exp((Ei-/i,)/T) 
1 + exp{Ei + qo - fij) _ 
'iT^'jiqo + mi — rrij — q^ /2mj 
2miq'^ 



In 



(8) 

(9) 
(10) 

(11) 



where and rrii are the chemical potential and the mass of baryons. In Eqs. ©-(IZI), 
Ee = qo + Ey, Ef.+ = qo + Ep, and E^. = —go ~Mp (see Fig. 1). Equations ®-(IZj) are valid 
for nonrelativistic and noninteracting baryons [22] . Below the nuclear density, this is a good 
approximation. It should be emphasized that Eqs. (jH))- (fTT|) differ from the corresponding 
equations in Reddy et al. in which the mass difference between nucleons is neglected. 
However, in order to investigate the rate of neutron decay, it is necessary to keep the mass 
difference in Eqs. ©-(IIII). 

n 

Analogous to the analysis in Reddy et al. |22], it is easy to obtain the previous results 
in the nondegenerate and degenerate limits of baryons. For instance, the results in the 



nondegenerate limit corresponding to Eqs. ©-((Zj) are shown below, 

1 f°° 

Xe-p ^ TT^m^nj, / dE,E,p,{E,-Qff,, (12) 

^TT Jq 



1 1'°° 

— \M\^n^ dE,E,p,iE, + QfU, (13) 



JTte 



27ra 

X(l-/e), (14) 

where Q = m„ — mp is the mass difference between neutron and proton, rij = 
2(mjT/2'7r)^/^ exp(r7j) is the number density of neutrons and protons in the nondegenerate 
limit, and rji = {fii —mi)/T is the reduced chemical potential. The above approximate rates 
are frequently cited in the literature to discuss the kinetic equilibrium for the /5-processes 
and the emissivity of neutrino emission, even though its validity should be checked carefully 

At me, iieiQ < T, the above reaction rates can be simplified further into 

X exp(?7p + /ie/T - Q/T) 

x[/(4) + 2/(3)g/r], (15) 
X exp{r]n - fJ^e/T) 

x[/(4) + 2/(3)g/T], (16) 
A„ ~ 1.63 X 2-'\-^-^\M\^ml;'-T^-'ml, (17) 

here /(3) = 77rVl20, and /(4) = 45C(5)/2, C(5) = 1.037 is the Riemann (-function. The 
definition oiUn) is I{n) = J^{e^ + l)^^x"'dx. This integral can be found in Gradshteyn 
and Ryzhik [24 1. 



In the elastic limit, Up in Eq. (fT^ is replaced by 

Vpn = inp-nr,)/il-e^"-^/^). (18) 

It is generally believed that electron capture rate under the elastic approximation in some 
sense introduces the effects of the degeneracy of baryons, so it should be more accurate 
than that under the nondegenerate approximation. Suppose that the nuclei are dissolved 
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FIG. 2: Electron capture rates on protons in the dense subnuclear matter out of the chemical 
equilibrium as a function of the baryon number density ne at different temperatures. The electron 
fraction is taken to be 0.3. The solid curves show the results from Eq. the dashed curves 
are the elastic results from Eq. (|18|) . the dot-dashed curves are the nondegenerate results from 
Eq. (HI. 



completely into nucleons at high temperature. Figure |21 shows the differences between our 
electron capture rate in the dense subnuclear matter out of the chemical equilibrium and 
the previous approximate results. The electron fraction Ye = (^e- — ne+)/n-B is assumed 
to be 0.3. The number density of particles at any degeneracy is expressed in terms of the 
Fermi-Dirac functions [25 1. The multidimensional integrations and the Fermi-Dirac functions 
are calculated using the mixture of Gauss-Legendre and Gauss-Laguerre quadratures {2^. 
From Fig. |21 it is evident that there are great differences between the elastic results and the 
general results (Eq. when baryons become degenerate. As shown in Eq. (fTHjl . the capture 
rates under the elastic approximation always decrease exponentially when nucleons become 
degenerate, which is qualitatively correct in the dense nuclear matter near the /3-equilibrium. 
However, this conclusion is not correct obviously in the dense subnuclear matter out of the 
chemical equilibrium. The elastic approximation is also good when fip ^ yU„ because the 
favored energy transfer is zero. In short, the elastic approximation underestimates the 
electron capture rate, when baryons become degenerate. 
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III. CONDITIONS FOR KINETIC EQUILIBRIUM. 



A. Chemical equilibrium condition for neutrino trapping. 

If neutrinos are trapped in a system, as what happens in the interior of a PNS, the inverse 
reactions corresponding to beta processes (H}-© by absorption of neutrinos can take place: 

e~ +p ^ n + Ue, (19) 
e+ + n ^ p + i^e, (20) 
n ^ p + e~ + i^e, (21) 

and their rates are equal to these of the corresponding beta reactions when the system 
reaches the chemical equilibrium. Because both the photons and neutrinos are trapped, the 
chemical equilibriums 7 + 7 ^ e"*" + e~ ^ z/g + t'e give(/i-y = 0), 

fie+ = -/^e-, (22) 
f^p, = -f^u,- (23) 

Based on the standard arguments in the theory of thermodynamics, the chemical equilibrium 
of any one of the above reversible reactions (fT^ - ()2H) gives the same chemical equilibrium 
condition for both cold npe~ and hot npe^ gases, 

fJ^n + fJ'U, = fip + fie- (24) 

If the chemical potential of the trapped neutrinos is zero, the above equation reduces to 

A^n = yUp + fie, (25) 

and the corresponding number density of the trapped neutrinos is written as 

T3 p x'dx 3T3 

where ({3) ~ 1.202. It is clearly shown in Eq. ()26p that the number density of the trapped 
neutrino is not zero at all when their chemical potential is zero. Only when T — 0, n^^ = 
rty^ — * 0. Therefore, Eq. (f^H|l can also be understood as the chemical equilibrium condition 
for cold npe~ gas under /3-equilibrium when neutrinos can escape freely from the system. 
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B. Kinetic equilibrium condition for neutrino escaping. 



In principle, if the matter is transparent to neutrinos which carry energy away, the equi- 
hbrium of the reactions CQ)-® can not be treated as a chemical equilibrium problem 20|. 
Suppose that the dynamic time scale of the system under consideration is greater than that 
of the reactions (P)-©, the general condition for the /3-equihbrium is given by 

= ^e+n{l^n, f^p, i^'ei T) + Xnif^'n, f^p, l^ei T). (27) 

Given n-Q and T, there are two additional trivial conditions. One is the conservation of the 
baryon number, 

n„(/in, T) + np(/ip, T) = ne, (28) 
the other is the charge neutrality, 

ne-(/ie,T) - ne+{ne,T) = np{^p,T). (29) 

Consequently, the chemical potentials of neutron, proton and electron {fin, fJ'p, fJ^e) can be 
determined by a set of closed equations, i.e. , Eqs. 



IV. ANALYTIC /3-EQUILIBRIUM CONDITIONS FOR NEUTRINO ESCAPING. 
A. /3-equilibrium condition for cold ripe gas. 

For ideal, cold npe~ gas, we can set the chemical potential of the neutrinos in Eq. (j24|) 
to zero because their number density is zero at T = (see Eq. (|26j) ). therefore, the chemical 
equilibrium conditions for the ideal, cold npe~ gas is Eq. (j2Sl), i-e- , fJ'n = fJ'p + IJ-e 13 • 

In the following, we will re-derive the chemical equilibrium condition for the cold npe~ 
gas Eq. from the viewpoint of the equality of the reaction rates of the beta processes. In 
the interior of an old neutron star, the electrons are degenerate {Ee ^ T), thus the number 
density of positrons are exponentially depressed. As the typical energy of emitted neutrinos 
is of order the temperature, the energy and the momentum of the emitted neutrinos can be 
neglected comparing to the Fermi energy of particles, i.e. , Ep + = En, and Pp+pl = Pn- 
So the kinetic equilibrium requires Ag-p = A„, therefore. 



2 n 



1=1 



(2vr) = 



(27r)^|M|25(4)(P, + Pp-P„; 



[/p(l-/„)/e-(l-/p)/n(l-/e)] 



2 n 



i=l 



(27rr|Mp5W(Pe + P,-P„, 



(30) 



In the above derivations, we use the even property of the Dirac-delta function, i. e. , 5*^^-* (Pg + 
Pp — Pn) = ^''^''{Pn — Pe — Pp)- It is evident that Eq. (IHUj) eventually results in the well-known 
chemical equilibrium condition Eq. ()25|) which is generally used to determine the equation 
of state of the dense matter in the interior of an old neutron star. 



B. Analytic steady state equilibrium condition for hot npe^ gas. 

It is well known that in a system with plenty of positrons, the pair must not be 
degenerate (T > i?e), if not, the number density of positrons will decrease exponentially, 
therefore, exp((i?e — f^e)/T) > 1. Before we derive the analytic dynamical equilibrium 
condition for such system from which the emitted neutrinos can escape freely, we first explain 
why the condition is not the well-known result of Eq. ()25|) . Suppose that Eq. is satisfied, 
for the reason described in §111. A., the chemical potential of the trapped neutrinos is zero, 
thus the ratio of the number density of the trapped neutrinos to that of the electrons is about 
T^/E^ > 1, if the electrons are not degenerate! This conclusion is obviously contradictory 
to the precondition that the system is transparent to neutrinos. 

In such system, comparing to the rate of positron capture by neutrons, the rate of neutron 
decay could be neglected before reaching the degenerate limit. Thus the /3-equilibrium 
condition for hot npe^ gas is 

m 

^e-p = K+n^ (31) 

if the hot electron-positron plasma with nucleon is transparent to neutrinos. Another good 
approximation is the elastic approximation, that is, the energy of the emitted neutrinos is 
of order that of the captured e^, i.e. , E^, ~ Ee, Ep ~ Ee+, and thus En — Ep. 
Under these approximations, the kinetic equilibrium condition Eq. (jHTj) gives 

= Ae-p — Ag+n 
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2/n 



i=l 



(27r)^|M|V^)(Pe + Pp-Pn-P., 



[/p/e(l-/„)-(l-/p)/e+/«] 
- fn)fe+ 

g(i?e-A«e)/T _^ X 



(27r)^|Mp5W(Pe + Pp-Pn-P. 



,{(£^p-E„)-{/ip-/i„))/T 



2/n 



j=i 



(2vr)= 



(27r)^|Mp5W(Pe + Pp-Pn-P.. 



(32) 



Therefore, from Eq. (fH^ we obtain the beta-equihbrium condition for the plasma with 
nucleons 

/Un = /ip + 2/ie. (33) 

It should be emphasized that during the above derivation, it is not assumed whether the 
baryons are degenerate or not. On the other hand, Eq. is still valid under the degeneracy 
of baryons, which is neglected completely in the approximate reaction rates at the beginning. 
However, it should be pointed out that wherever baryons are degenerate in astrophysics, 
neutrinos are generally trapped and likely thermalized. In these post-neutrino trapping 
environments, the analytic steady state equilibrium condition we derived here is not valid, 
because the effects of neutrino trapping are not included in our consideration. It is not 
a surprise to notice that the analytic equilibrium condition Eq. can also be drawn 
in the completely nondegenerate limit. Setting the equality of Eq. (fT3|) and Eq. (fTBj) . the 
equilibrium condition is also obtained. In the following, we have another heuristic but not 
very strict derivation. If all particles are nondegenerate, we have 



\e-p oc rie-Up oc exp{r]p + r]e), 
Xe+n OC ne+n„ OC exp(?7„ - ?7e). 



(34) 
(35) 



Setting the equality of Eq. and Eq. (jH^ . the analytic condition Eq. (jHHI) is obtained 
again. 

The advantage of having the analytic equilibrium condition at hand is obvious. For 
instance, we can derive some useful formula in the nondegenerate limit of baryons. In such 
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limit, using the Saha equation gives the ratio of neutron baryon number density to that of 
! I 

proton {e.g. [26]). 

— = 1 - — = e T =eT. (36) 

At the last step, the analytic equilibrium condition Eq. (jH!^ is applied. On the other hand, 
for the relativistic e^, there exits an exact expression for the electric charge density in terms 
of the chemical potential of electrons {e.g. |27|). 



T3 



(37) 



37r2 

Substituting Eq. into Eq. (jHTj), we have only one simple equation to determine the 
electron chemical potential /ig. 
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37r2 



0. (38) 



1 + g(2Ate-Q)/T 



After that, the electron fraction Y^. can be directly calculated from Eq. (j37p . 
At iie/T < 1 and (2/ie — Q)/T < 1, Eq. (jHH|) can be simplified further, 

_ 1 (i + o.sg/T) 

2(l + 1.5nB/T3)- ^^^^ 

Equating the rates of e~ and e"*" captures in the region of mild degeneracy, Beloborodov 
obtained a similar result to Eq. (j39|) i see Eq. (11) in his paper). In a word, using the well 
known results and our new result Eq. ()33p. we re-derive a similar previous result in a different 
way. From this consistence, it turns out that Eq. ()33|) is correct and the approximate formula 
Eq. ()38|) is more accurate than the previous result. 

In the following, we check the validity of Eq. ()33|) by numerical calculations. As before, 
suppose that the nuclei are dissolved completely into nucleons at high temperature. Given 
n-Q and T, the electron fraction Yf. can be determined by two different sets of the equilibrium 
conditions. One set consists of Eqs. (fTTfl - lj^ . In Eq. (j2Zj), the [3 reaction rates are calculated 
based on the reaction rates Eqs. ©-©j or the approximate rates Eqs. (fT ^ - (fT^ . respectively. 
The other set consists of our analytic equation Eq. (jHHj) . and Eqs. (f^ - (f^ . Figure El shows 
Ye and the reduced chemical potentials rji {i = n, p, e) versus the baryon number density at 
different temperatures. In the upper panel of Fig. Q, the solid lines show the results from 
Eqs. (j27j)-(j2ni) and Eqs. (jSI)-© for the P reaction rates. The dot-dashed curves are from 
Eqs. f07 |) - (j^ . but the approximate reaction rates in the nondegenerate limit Eqs. (fT^ -(fT ^ 
are used. The dotted lines represent the results from Eq. (fH^ and Eqs. (| ^ -(| ^ . It is 
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0.1 

FIG. 3: Upper panel: Electron fraction Yg under /3-equilibrium versus the baryon number density 
n-Q at different temperatures. The soHd curves show the results from Eqs. ©-©) the dotted 
curves are from the analytic /^-equilibrium condition Eq. and the dot-dashed curves are 

from the nondegenerate approximation of baryons. Lower panel: The reduced chemical potentials 
of electrons, neutrons, and protons rji as a function of the baryon number density at different 
temperatures. From top to bottom, the lines correspond to the results of T = 6 x 10^, 10^'', 
5 X 10^°, 10"K. 
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FIG. 4: The neutrino emissivity versus the baryon number density at different temperatures. 
The lines are the same as those of the upper panel in Fig. |31 

evident that the resuhs from the analytic condition are almost consistent with the general 
numerical results. As expected, the difference occurs only in the regime where electrons 
become degenerate, beyond which, the condition for the analytic formula Eq. ()33|) is not 
satisfied any longer and the contribution from the neutron decay can not be neglected. 
As shown in the lower panel of Fig. neutrons become degenerate as electrons do. The 
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particles become degenerate if their reduced chemical potential exceeds the temperature. 
For electrons and nucleons, their degenerate number densities are estimated to be 
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(ji \ 3/2 
I^ttkJ f^"'- (41) 

In the same regime, there are great differences (of several orders) between our results and 
the previous results in which the degeneracy of nucleons is completely neglected. In any 
case, Fig. IHl evidently shows that the results from the analytic condition are much more 
accurate than those from the approximate rates in all the parameter regions. 

The total neutrino emissivity under /3-equilibrium is shown in Fig. |3] Compared with the 
general numerical results, both approximate methods overestimate the rate of the neutrino 
emission because neglect of the degeneracy of particles increases the phase space for the 
relevant reactions. It is clearly shown in Fig. |3 that the effect of the neutron degeneracy is 
much more important than that of electrons, if such conditions are satisfied. 



V. CONCLUSIONS AND DISCUSSIONS 

n 

In this work, using the structure function formalism developed by Reddy et al. [22|, we 
derive the rates of the (3 processes which including captures and neutron decay in the 
dense subnuclear matter for all participating particles are at arbitrary degeneracy. For this 
purpose, the difference between the mass of neutron and proton is kept in the structure 
function (see Eq. (|TT|l ). which is neglected originally in Reddy et al. [22|. Comparing to 
our reaction rates, the previous approximations have been checked. Especially, the electron 
capture rate under the elastic approximation which is put forward to include the effects 
of the degeneracy of baryons differs dramatically from our results when baryons become 
degenerate and before the dense matter reaches the /3-equilibrium. Generally speaking, 
the elastic approximation underestimates the electron capture rate when baryons become 
degenerate. 

Most of interest, we find an analytic steady state equilibrium condition for plasma 
with nucleons, if the system is transparent to the emitted neutrinos. This result is valid 

14 



when electrons are nondegenerate or mildly degenerate, but it has nothing to do with the 
degeneracy of baryons. Basically, it is a good result for plasma with nucleons, because 
when electron becomes degenerate, the number density of positron decreases exponentially, 
then the system can not be called as plasma. If the nucleons are nondegenerate, we 
further obtain only one simple equation for determining the electron chemical potential and 
electron fraction, which shows the great advantage of our analytic equilibrium condition. 

So far, there are four analytic beta equilibrium or steady state equilibrium conditions 
for npe^ gas in different astrophysical circumstances. They are summarized as follows: 
(1), If the neutrinos are completely trapped in a system, such as in the interior of a PNS, 
/^n + /^z^e = /^p + /^e- (2), In the cold npe^ gas under beta equilibrium, /i„ = /Up + /Xg- (3), If 
the neutrinos are partially trapped, and the chemical potential of the trapped neutrinos is 
zero {^i,^ — 0), /In — IJ'p + IJ'e- (4), If the neutrinos can escape freely from the system with 
plenty of e^ pairs, //„ = //p + 2/ie, which is the main result of this work. 
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